Indonesian Journal of Electrical Engineering and Computer Science 
Vol. 17, No. 2, February 2020, pp. 662~670 
ISSN: 2502-4752, DOI: 10.1159 1/ijeecs.v17.i2.pp662-670 Oo 662 


Multi cases optimal reactive power dispatch using evolutionary 
programming 


Rahmatul Hidayah Salimin’, Ismail Musirin’, Zulkiffli Abdul Hamid’, Afdallyna Fathiyah Harun’, 
Saiful Izwan Suliman*, Hadi Suyono’, Rini Hasanah’ 
'234>Raculty of Electrical Engineering, Universiti Teknologi MARA, Malaysia 
° "Department of Electrical Engineering, Brawijaya University, Indonesia 











Article Info ABSTRACT 

Article history: Evolutionary Programming (EP) is one of many types in Evolutionary 
; Computation (EC) that used for optimization process. EP technique is used to 

Received Jun 10, 2019 find the optimal reactive power dispatch (ORPD) since it is one of the 

Revised Aug 13, 2019 accessible options schemes that can be used on the system as a reactive 

Accepted Aug 27, 2019 power support. Sometimes, it is not necessary to operate all generators in 


order to perform ORPD to in achieve the objectives. Also, increment of 





reactive power load to the system will cause voltage decomposes with the 
Keywords: increase in transmission loss in the system. Therefore, the proposed method 
decides the best grouping of generators that should be operated in system by 
bearing in mind the transmission loss reduction. ORPD will be used to 
. minimize the transmission loss as well the increasing reactive power loading. 
Multi cases ORPD This method conducted on IEEE 30-bus Test System with multi cases 
scenario. The best combination of operating generators determined and the 
transmission loss after optimization is smaller compared to the transmission 
loss before optimization resulted. 


Evolutionary programming 
Loss minimization 


Copyright © 2020 Institute of Advanced Engineering and Science. 
All rights reserved. 


Corresponding Author: 


Rahmatul Hidayah Salimin, 

Faculty of Electrical Engineering, 
Universiti Teknologi MARA, 

40450 Shah Alam, Selangor, Malaysia. 
Email: rhidayah @uitm.edu.my 








1. INTRODUCTION 

Power system is a composite network utilized for creating and transmitting electric power [1]. 
It needs a requirement to maintain the voltage stability in order to make the operation smoothly run. In power 
system operation, the power loading is not always fixed and usually it is increasing according to the demand. 
Along the increment of load demand, the voltage will slowly decrease with the increasing transmission loss 
in the power system. The greater part of the load demand included reactive power impacts because of 
changes in burden like motor installing, transformers and generator blackouts as opposed to the real power 
effect [2]. It can lead to insecure power system operation condition since the power system is required to 
work with utilization of minimal resources giving greatest security and reliability [1]. In order to avoid this 
problem, a solution should be considered to sustain the increment of reactive power loading and thus the 
transmission loss will be minimized as well the maximized voltage deviation. Towards settling this problem, 
it is a need to manage the reactive power resources in power which can be solved using the ORPD. 

ORPD has a role that very significant in power system in order to ensure the economy, security and 
stability of a power system in secure condition. The main objective of the ORPD is to minimize the system 
real power loss [3]. Even ORPD is well known for transmission loss minimization, it also can used in 
appropriate scheduling of reactive power and its dispatch in the system [4]. Scheduling of reactive generators 
by determining the best combination operating system is suitable for decision making tool in order to use 
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certain generator in minimising the transmission losses and maximising the voltage deviation. ORPD is 
originally a sub-issue of optimal power flow (OPF) estimations. It has been defined as a non-linear, 
multi-uncertainly, multi-constraint, multi-minimum and multi-objective optimization [5]. OPF was first 
introduced by Dommel and Tinney and developed by Carpentier in 1970s [6]. ORPD is a reactive power 
support through existing generator bus, synchronous condenser, capacitors, static compensators and tap 
charging transformer. This problem denotes optimal settings of control variables to minimize a certain object 
while satisfying equality and inequality constraints [7]. 

Numerous classical methods with excellent convergence characteristic have been developed and 
applied to solve ORPD problem in the past decades, such as Interior Point Method (IPM), Linear 
Programming (LP) and Quadratic Programming (QP). These methods also can be as conventional 
optimization techniques or mathematical optimization techniques. However, there exist some limitations and 
challenges for these methods to handle non-differentiable objectives functions and discrete variables in 
ORPD. LP method requires that the objectives function and constraints have linear relationship, which may 
lead to loss of accuracy [1]. Besides that, the Gradient and Newton methods which also the methods of 
mathematical optimization techniques suffer from the difficulty in handling inequality constraints [1, 8, 9]. 
Zhu and Xiaong proposed a new approach to study the ORPD using a Modified Interior Point (MIP) method 
to minimize the system power losses and to penalize any new ORPD utilization while Granade et al. 
approached a decentralized approach based on Lagragian decomposition method for solving ORPD problem 
in multi-area power systems however still these classical methods suffer from many drawbacks [10, 11]. 
Therefore, these numerical-based techniques are not appropriate for taking care of the ORPD issues. 

Hence, more modern optimizations methods were proposed in order to improve the drawbacks of 
the classical methods in finding the ORPD solutions. Aruna and Devaraj has proposed an improved genetic 
algorithm which permits the control variables to be represented in their natural form have found the result 
that this method could minimize the losses and suited for solving the mixed integer optimization problem [8]. 
JAYA algorithm was proved to be well consistent, robust, stable and efficient in determining the optimal 
location of the control variables for minimizing the active power loss for all the test cases, while satisfying all 
the constraints simultaneously using ORPD [12]. Khazali and Kalantar have proposed Harmony Search 
Algorithm (HSA) for determination of the global or near global optimum solution. It has the advantage in 
decreasing transmission loss, voltage deviation and increasing voltage stability margin compare to simple 
genetic algorithm (SGA) and particle swarm optimization (PSO) [7]. Mukherjee have proposed chaotic krill 
herd algorithm (CKHA) to minimize the losses and minimize total voltage deviation. It is a novel metal- 
heuristic search approach which is introduced by Gandomi and Alavi in 2012 based on the simulation of the 
herding behavior of krill swarms in nature [13]. Besides that, exchange market algorithm (EMA) based on 
the behavior of shareholders in stock market also used to solve ORPD issues [14]. There are many more 
algorithms such as Sepiida algorithm (SA) [15], Modified Monkey Optimization (MMO) algorithm [16], 
Advanced Fireworks Algorithm (AFA) [17], Multi-Verse Optimizer (MVO) [18], Crow Search Algorithm 
(CSA) [19], Wolf Search Algorithm(WSA) [20], Rapid PSO (RPSO) [21], Gray Wolf Optimizer 
(GWO) [22,] and many more have been proposed in order to achieve the objectives of the ORPD besides the 
methods that have been mentioned above. 

This paper proposed Evolutionary Programming (EP) algorithm for solving ORPD issues. This 
method is used to determine the reactive power required to be dispatched by generators in the system and to 
determine the best grouping of generator that should be operated. The objective function of this tactic is to 
minimize the transmission loss along the increasing of power loading. EP optimization engine has been 
established to apply the reactive power dispatch processes. This method was verified on IEEE 30 bus system. 


2. METHODOLOGY 
2.1. Evolutionary Programming Technique 

Evolutionary Programming (EP) is a stochastic optimization technique which is inspired by the 
theory of evolution by means of natural selection and was developed by Lawrence J. Fogel in in 1960 for 
prediction of finite state machine [23, 24]. There are a few types of EP such as classical EP, meta EP, 
and adaptive EP which every type of EP have different mutation scale. EP finds for ideal solution by 
developing a population of candidates or parent population over a number of generations. Another population 
also known as offspring is shaped from the current population (parent) by executing the mutation scale 
during every generation. The mutation operator delivers afresh solution by perturbing every segment of the 
present solution by a random amount. The fitness of every candidate solution dictates its strength. The fitness 
is evaluated from objective function of the optimization problem. In the selection process, every individual 
from same population will for new generation. The new generation will select the individual that gained the 
most numbers of successes. The individual that obtained the least numbers of win will be eliminated. These 
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processes involved in EP will lead to achieve the global optimal solution. The flowchart of EP is shown in 
Figure 1. 
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Figure |. Flowchart of EP operation 


2.1.1 Random Number Generation 

EP operates with a set of population that contains the number of individuals or also known as parent 
population. At first, the parent population can be produced by generating random number using a uniform 
distribution number, where: 


ie A AAG Ab ce AG | () 


i= 1,2,3,4,5...n, where n is the population size with values between the control variable range which is from 
Amin tO Amax In the optimisation process. These random numbers are the new value of reactive power to be 
dispatched by the generators. 


2.1.2 Assign Load Bus for Test and Run Load Flow 

The number of variables, x from (1) represents the number of generators will be used. The symbol 
of x is unfixed. The reactive power on the chosen loads will be represented by parameter in estimating the 
maximum loading. 
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2.1.3 Set Constraint 

A few constraints are set toward the start to enable EP to create just random numbers that fulfil 
some predetermined constraints. If constraint is satisfied, the random number will fill the population. If not, 
new random number will be created. 


2.1.4 Fitness 1 Calculation 

The objective function will be considered and taken as fitness equation. For this study, transmission 
losses and minimum voltage are objective function which needs to be minimized and maximized within the 
range by conducting the load flow program. Fitness | calculation is the values of the objective function 
calculated using the first generated population known as parent. 


2.1.5 Mutation 

Mutation is a process to produce the offspring by some random modification of individuals in parent 
population. This process is implemented based on the (2). The mutation scale is manually adjustable in order 
to achieve better convergence. Small value of mutation scale causes small search step which leads to fast 
convergence of EP and vice versa. 





Xiemj = %Xiy + N (0.20 shits. Sita) (44) (2) 
Where: 
Xi+mj = mutated parents (offspring) 
Xij = parents 
N = Gaussian random variable with mean and variance 
Bb = mutation scale, 0 to | 
Xjmax = Maximum random number for every variable 
Xjmin = Minimum random number for every variable 
fi = fitness for random number 
fax = Maximum fitness 


2.1.6 Fitness 2 calculation 
Fitness 2 calculation is a recalculation process of objective functions using the mutated parents or 
offspring. 


2.1.7 Combination 
In this process, all individuals from parents and offspring will be compiled into new large 
population. 


2.1.8 Selection 

This process is used to determine the survivors for next generation where the survivor with better 
fitness values in the population have big chance of contributing offspring for next generation. Usually, the 
first half of new large population will be retained as the new individuals while the remaining will be 
discarded. 


2.1.9 Convergence Test 

In order to decide the stopping point of the generation process, this step needs to be performed. For 
decision making, the values of maximum fitness and the values of minimum fitness need to be compared and 
the difference should be less or equal to 0.0001. The program will stop if the solution converges and if not, 
Fitness | will be calculated again. 


2.2. Optimal Reactive Power Dispatch (ORPD) using EP 

The total lost minimisation was assigned as the objective function and used to implement the ORPD 
in the system. The variables of this study are the reactive powers that injected on the generator buses for 
decreasing the total transmission loss of the test system. EP was used to decide the ideal reactive power to be 
transmitted by the generator bus that used. MATLAB was used in developing the EP program. The ORPD is 
applied in order to decide the best combination of generator that should utilized as the responsive power 
compensator backing to the framework and thus improving the power framework security by considering the 
transmission loss. This can be a decision making guideline that can help the power system operator which is 
required by service organizations. 
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Random numbers that generated during initialization process will be defined as the injected reactive 
power of the generator buses in the system. The number of generator bus used in the system would influence 
the number of variables excluding the slack bus. The reactive power at generators 2, 5, 8, 11, and 13 would 
be controlled and were represented by five variables by considering the objective; to minimize power loss 
with load subjected to the system. These five variables are relegated as the reactive load with negative sign 
showing the reactive power are really infused at the specific buses. A constraint is set as the loss should be 
less or equal to loss_set. Loss_set here is a loss before the optimization started. 


2.3. Test System 

ORPD using EP was tested on the IEEE 30-bus Reliability Test System for optimising the total loss 
in the system. There are 24 load buses and 6 generator buses with 41 interconnected lines in this system. 
Multi cases scenario are created with various combination of five generator buses excluding the slack bus 
which used as the control variables for the ORPD. Selected generators are indicated with the tick sign. 


3. RESULTS AND ANALYSIS 

The study was verified at an assortment of combinations of generators and at various value of power 
loading of bus 15 on the IEEE 30-bus RTS to distinguish the outcome to the transmission loss and the 
minimum voltage. There was one constraint assigned in every grouping of producer before the RPD is 
executed with EP. The optimization operation is implemented at the various loading of bus 15 from 20VAR 
until 14OMVAR with increment of 20MVAR. The maximum loading of bus 15 is ISOMVAR. This loading 
bus was selected randomly. The result for transmission loss as the objective function is tabulated in the tables 
as following. The test system parameters values for pre-RDM are set such as mutation scale, 0.0015 and the 
population with 20 individuals. 


3.1. Multi Cases Generator Combination ORPD 

31 possible grouping for performing the ORPD procedure is applied to 31 possible grouping with 
five variables were used. The generator 2, 5, 8, 11, and 13 are utilized as control variables to dispatch the 
reactive power into the power system in order to observe that the minimum transmission loss is obtained. 
These tables show the ranking of generator grouping based on the transmission loss. The outcome of 
transmission loss minimisation when the loading bus 15 is varies between 80Mvar to 140Mvar is shown in 
Table 1-4. 

For 80Mvar loading bus, as shown in Table 1, the result shows that when only generator 5 works in 
transmitting the reactive power into the power framework, the obtained transmission loss is the most. The 
transmission loss is only decreased from 27.8308MW to 27.7920MW which is 0.12% reduction from 
transmission loss before optimization. At the same moment, the voltage profile is increased from 0.6709p.u 
to 0.6882p.u. On the other hand, the minimum broadcast loss is gained when generator 2, 8, 11, and 13 are 
forwarding the reactive power into the power framework. It changes from 27.8308MW to 24.7412MW with 
losses as much as 3.0891MW and 11.29% of reduction from transmission loss before optimization while the 
minimum voltage has increased from 0.6709p.u to 0.7724p.u. 

Table 2 shows the result of minimised transmission loss when the loading of bus 15 is IOOMVAR. 
From observation, the biggest transmission loss is gained when producer 5 delivers the reactive power into 
the power system. It changes from 32.6047MW to 32.3583MW with a difference of 0.2464MW. It has the 
highest percentage of transmission loss reduction from before optimization with the value of 14.87%. 
Simultaneously, the minimum voltage has slightly improved from 0.6422p.u to 0.6556p.u. However, the 
slightest transmission loss is gained when generator 2, 8, 11, and 13 are operating in delivering the reactive 
power into the power framework. The transmission loss is moderately decreased from 32.6047MW to 
27.7548MW with a loss of 4.8499MW with the voltage increased from 0.6422p.u to 0.7591p.u. It has 0.76% 
of transmission loss reduction from before optimization. 

Result that tabulated in Table 3 is about the transmission loss minimisation when the loading bus in 
100MVAR. During generator 2, 5, 8, 11, and 13 operating, the obtained transmission loss is the least one 
while the highest one is during generator 5 operates to transmit the reactive power into the power system. For 
the least transmission loss, it changes from 39.5358MW to 31.9114MW with the voltage increased from 
0.6154p.u to 0.7754p.u while for the most transmission loss is from 39.5358MW to 39.3988MW with the 
voltage incremented from 0.6154p.u to 0.6387p.u after performed optimization. The least transmission loss 
has the most percentage of the transmission loss reduction from before optimization with the value of 19.28% 
while in contrast, 0.35% for the most transmission loss. 
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Table 1. Transmission Loss at bus Loading 80Mvar 
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Table 3. Transmission Loss at Bus Loading 120Mvar 
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Table 4. Transmission Loss at Bus Loading 140Mvar 
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Table 4 shows the result of transmission loss minimisation when the loading of bus 15 is 
140MVAR. The highest percentage of reduction from transmission loss before optimization is when the 
generator 2, 5, 8, 11, and 13 operate in sending the responsive power into the power framework. 
The transmission loss is observed decreasing from 53.2168MW to 38.0275MW with the losses as much as 
15.1893MW. At meantime, the voltage profile incremented from 0.5686p.u to 0.7501p.u. On the other hand, 
the lowest percentage of reduction from transmission loss before optimization is when generator 2 does its 
work to dispatch the reactive power into the power system where the transmission has most loss. It is seen 
that the transmission loss is decremented from 53.2168MW to 52.3333MW with the least loss. While its 
voltage slightly increased from 0.5686p.u to 0.5821p.u. 


3.2. The Best 3 of Selected Operating Generator Combination 
The percentage of transmission loss is obtained by using the following equation: 


transmission loss before - transmission loss after 
—« eS _ ee— ll x 100 (3) 


transmission loss before 


Table 5 shows the best three combination of operating generator in minimising the transmission loss 
at different loading value of bus 15. These three combinations are selected from 31 possible combinations of 
generator 2, 5, 8, 11, and 13. From this table, it is observed that when the loading is increasing, the 
percentage of transmission loss reduction by using RPD is also increasing. It also shows that it is not 
necessary to operate all generators for minimising the transmission loss. For example, only three generators 
are used to achieve the better minimised transmission loss at loading of 20MVAR. It showed that by 
determining the combination of which generators is important in giving the better result for achieving the 
objective. 


Table 5. The Best Three Combination of Operating Generators 
Loading Percentage | Optimal reactive Power Dispatch (MVAR) 
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4. CONCLUSION 

ORPD for transmission loss minimisation using EP as the optimisation approach was presented. 
Transmission loss is decided to be the fitness function and must be minimised along the increment of power 
loading at bus 15. The result obtained showed the objective to minimise the transmission loss is achieved. At 
the same time, the best three combination of operating generators was determined based on ranking of the 
transmission loss. This result from the multi cases combination of operating generators can be used as a guide 
in helping the system operator in scheduling the generators and can avoid all generators to operate at one 
time. Increment of the loading bus showed the increasing of transmission loss. However, this issue can be 
overcome by ORPD using EP. In the conclusion, ORPD using EP is very efficient and can give a very good 
performance in minimising the transmission loss. 
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